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ABSTRACT 

Wc have performed hydro-dynamic and radiative transfer calculations of a photoevap- 
orating disc around a Herbig Ae/Be star to determine the evolution and observational 
impact of dust entrained in the wind. We find that the wind selectively entrains grains 
of different sizes at different radii resulting in a dust population that varies spatially 
and increases with height above the disc at radii > 10 AU. This variable grain popula- 
tion results in a 'wingnut' morphology to the dust density distribution. We calculate 
images of this dust distribution at NIR wavelengths that also show a wingnut mor- 
phology at all wavelengths considered. We have also considered the contribution that 
crystalline dust grains will have in the wind and show that a photoevaporative wind 
can result in a significant crystallinity fraction at all radii, when the disc is edge-on. 
However, when the disc's photosphere is unobscured, a photoevaporative wind makes 
no contribution to the observable crystallinity fraction in the disc. Finally, we conclude 
that the analysis of extended emission around edge-on discs could provide a new and 
independent method of testing photoevaporation models. 

Key words: accretion, accretion discs - circumstellar matter - planetary systems: pro- 
toplanetary discs - stars: pre-main-sequence. 



1 INTRODUCTION 

Protoplanetary discs are ubiquitous around young stars, 
an inevitable consequence of angular momentum conserva- 
tion during star formation. Young stars accrete their mass 
through discs which then provide the reservoir of material 
from which planets may form. Any understanding of the pro- 
cesses involved in star and planet formation requires knowl- 
edge of how discs evolve and are finally dispersed. Obser- 
vations of disc evolution have typically focused on thermal 
emission from the dust grains in the disc through the analy- 
sis of their continuum spectral energy distribution (SED) at 
several different wavelengths. While empirical analysis and 
numerical modelling of the SED of discs around young stars 
have proved extremely useful to study this important phase 
of the star and planet formation process, a number of well 
known degeneracies conspire to render the interpretation of 
the data non-unique. To remedy these shortcomings recent 
observational efforts are making use of various techniques to 
obtain spatially resolved information of many sources. 

Herbig Ae/Be stars are more luminous than their lower 
mass T-Tauri counterparts and provide an opportunity to 
observe large scale morphology in resolved images. Discs ob- 
served at nearly edge-on inclinations are particularly useful 
in studying the large scale structure since the star and inner 
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disc, which are much brighter than any extended emission, 
are obscured along the line of sight. Several edge-on discs 
around young stars exhibit extended emission above and 
below their mid-plane at NIR wavelengths (e.g. Padgett et 
al. 1999 & Perrin et al. 2006) and some have been further 
imaged through AO systems providing excellent spatial res- 
olution, as in the case of PDS 144N (Perrin et al. 2006). 

Photoevaporation (Hollenbach et al. 1994) is the mech- 
anism by which high energy photons (UV- X-rays) heat the 
surface layers of the discs to temperatures of order the escape 
temperature, forming a thermally driven wind from the disc 
out to large radius. Photoevaporation has been proposed 
as the dominant disc dispersal mechanism around low-mass 
stars (e.g. Clarke et al. 2001) and has been successful in 
explaining many of the observational statistics (Owen et 
al. 2010b). However there is only tentative evidence that 
the photoevaporative wind itself has been detected, through 
Nell and OI emission (Hartigan et al. 1995; Alexander et 
al. 2008; Pascucci et al 2009; Schisano et al 2010; Ercolano 
& Owen 2010) which probes the wind on a very local scale. 
In order to test photoevaporation, along with the various 
different heating mechanisms we need global observational 
diagnostics to validate the models. Edge-on discs provide 
the perfect opportunity to detect a disc wind on a global 
scale through light scattered by the dust grains that will 
inevitably be entrained by the wind. 

The origin of crystalline dust grains detected in proto- 
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planetary discs via infra-red (IR) spectroscopy (e.g. Bouw- 
man et al. 2001) and directly within our own solar sys- 
tem (e.g. Wooden et al. 1999) is still a matter of discus- 
sion. The source of all dust in star formation - the ISM 
- is inferred to be entirely amorphous (e.g. Kemper et al. 
2005). Spectroscopy of discs (van Boekel et al. 2005; Apai et 
al. 2005) indicates a non-negligible level of crystalline dust 
outside the crystalline radius, where grains are hot enough 
(T > 800K) to be thermally annealed and converted into 
crystalline grains. This discovery has lead to the develop- 
ment of disc models with radial mixing to allow crystalline 
grains that formed in the hot inner disc to be transported 
to larger radii (e.g. Morrill & Voelk 1984; Gail 2001, 2002; 
Wehrstedt & Gail 2002; Bockelee-Morvan et al. 2002; Dulle- 
mond et al. 2006; Hughes & Armitage 2010). Another sug- 
gestion put forward by Shu et al. (1996), is that crystalline 
grains could be carried outwards by a x-wind. Given that the 
crystallization radius of <1-2AU is smaller than the min- 
imum launch radius for a photoevaporative wind > 2AU 
around Herbig Ae/Be stars and lower mass stars, the pho- 
toevaporative wind cannot be the direct source of the crys- 
talline grains at larger radii, but could work in combination 
with radial mixing to produce the observed enhancement. 
Furthermore crystallinity around edge-on discs can be used 
to probe the source of the dust in the extended emission, 
since a disc origin will give rise to crystalline dust grains 
in the wind while an in-falling envelope will contain only 
amorphous grains. 

In this paper we present a simple model to investigate 
the imprint of photoevaporation on the extended emission 
observed around edge-on Herbig Ae/Be stars. We will com- 
pare our model to current observations of edge-on discs along 
with following the fate of crystalline grains in the wind as 
they are transported to large radius. In Section [2] we de- 
scribe the model, including the hydrodynamic and radiative 
transfer methods. In Section Owe present synthetic images 
obtained from the models. In Section [4] we describe the re- 
sults of crystallinity calculations for the wind. We compare 
our results to observations of edge-on discs in Section [5] and 
we summarise our main findings in Section [6] 



launching surface, it will be entrained permanently, allowing 
dust to be carried to very large distance from the star, since 
the drag force will dominate over gravity. 

We build a simplified model to test the possibility of 
extended emission from dust grains due to a photoevapo- 
rative wind, and we choose a set of assumptions that al- 
lows us to place an upper limit on this expected level of 
emission. Namely we assume that the entire dust popula- 
tion is able to reach the launching surface of the wind via 
some turbulent mechanism and ignore the effects of settling 
and grain growth. Our simple model ignores the details of 
the bound regions of the disc0 and thus ignores any emis- 
sion produced by this region, including only the wind itself 
and its contributions to emission. This simplification means 
our model will not be accurate near the mid-plane, where 
emission from dust in the disc's upper atmosphere domi- 
nates over the wind emission. Therefore, we are unable to 
reproduce the observed intensities and optical depths seen 
at several scale height above and below the mid-plane in an 
edge-on disc; however, here we are interested in the extended 
emission, which, in the absence of an infalling envelope is 
dominated by the photoevaporative wind. 

Our model naturally produces a 'wingnut' morphology 
similar to that seen in several scattered light images (e.g. 
Perrin et al. 2006), and produces a spatially variable dust 
distribution, due to the different maximum grain size that 
can be entrained along each streamline. The radial morphol- 
ogy of the streamlines then yields that, at a given cylindri- 
cal radius, the maximum grain size increases with height, 
resulting in a spatial variation of colour in the scattered 
light images. Furthermore crystalline grains entrained in the 
wind are transported outwards from regions of the disc with 
higher crystallinity fractions resulting in an enhancement 
of crystalline grains in the wind over the disc's underlying 
distrubtion. 

We can split the method for constructing scattered light 
images into three separate parts: (i) hydrodynamic calcula- 
tions of the photoevaporative wind; (ii) calculation of the 
dust profile distribution and crystallinity based on the hy- 
drodynamic solution; (iii) radiative transfer modeling of the 
dust distribution. 



2 MODEL 

Considering the force balance on a dust grain, it is easy to 
show that small dust particles are entrained in a photoe- 
vaporative wind, which carries them out to large distances. 
Takeuchi et al. (2005) showed that the drag force on a grain 
is approximately: 

m d p w vl 

rd (1) 

Pda 

where rrid, pd and a are, respectively; the mass, the density 
and the radius of a dust grain (assumed to be spherical) 
and p w and v m are the density and velocity of the wind. At 
large radius 0(z/.R > 1) the flow is approximately spherical 
implying p w v w falls off as 1/r 2 . Since gravity also falls off 
as 1/r 2 and v w increases monotonically with radius, then 
if a grain is still entrained several scale heights above the 



2.1 Model Set 

As the mass- loss rates are not well known for Herbig Ae/Be 
stars we have left the ionizing luminosity (which sets the 
mass-loss rates) for an EUV driven wind (Hollenbach et al. 
1994) as a free parameter and consider the effect of chang- 
ing mass-loss rates on the morphology, colour of the emis- 
sion and crystallinity distribution in the winds. We consider 
an EUV wind from a primordial disc around a 2.5M© star 
with a range in ionizing luminosity from 10 41 to 10 45 s^ 1 in 
steps of 1 dex. These luminosities correspond to mass-loss 
rates in the range 10~ 10 — lO~ 8 M0yr -1 and are similar to 
values calculated by Alexander et al. (2004). In order to cal- 
culate images we compute scattered light images at 1.6, 2.1 
& 3.8/im with band passes corresponding to the H, K' & 
L' bands. In order to consider the relative contributions of 



1 Throughout this work we use R, <f>, z to refer to cylindrical co- 
ordinates and r, 0, (p to refer to spherical polar co-ordinates 2 In the sense that T <C T cs 
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crystalline to amorphous dust we take the calculated crys- 
tallinity fractions in the disc of Dullemond et al. (2006) for 
an evolved Herbig Ae/Be star. We assume that this disc 
crystallinity fraction is the same at the base of our photo- 
evaporative wind, such that the crystallinity fraction along 
the streamline can be calculated by following the thermal 
evolution of the dust entrained on that streamline. 



2.2 Hydrodynamic photoevaporative wind 

A photoevaporative disc wind is a thermally driven hydro- 
dynamic wind occuring when the disc's surface is heated to 
temperatures of order the escape temperature, allowing it 
to launch a freely expanding wind. While the wind driv- 
ing source for lower mass (T-Tauri) stars is likely to be 
X-rays (Owen et al. 2010a, Ercolano & Owen 2010, Er- 
colano & Clarke 2010, Owen et al. 2010b) the wind driv- 
ing source around intermediate mass stars has not yet been 
thouroughly investigated. X-ray photoevaporation may still 
occur to some degree; however, the lower Lx / Lboi ratio of 
Herbig Ae/Be compared to T-Tauri stars and their higher 
EUV fluxes, may reduce the role of X-rays in driving the 
wind. We adopt here the EUV driven wind of Hollenbach et 
al (1994) and hydrodynamic solution of Font et al. (2004), 
which is a simple and scalable hydrodynamic solution, al- 
lowing us to consider a wide range of parameter space, 
something not possible with more complicated FUV mod- 
els (Gorti & Hollenbach 2009) or X-ray models (Owen et 
al. 2010a). The simplified EUV treatment is suitable for the 
purpose of this work which aims at being the first approach 
in studying the qualitative aspects of scattered light emis- 
sion from a disc wind. We have repeated the calculation 
of Font et al. (2004) and Alexander (2008) and we refer the 
reader to these papers for a detailed description of the model 
setup, while the basics are summarised below. 
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Figure 1. EUV wind structure: the colour-map shows the density, 
while the arrows show the direction and magnitude of the gas 
flow, note that just after launch the flow is very close to being 
spherical. 



where n g is the density at R g which was determined t hrough 
the numerical calculations of lHollenbach et alj (|l994l ) to be: 



2.8 x 10' 



( * 



\10 41 s 



1/2 



-3/2 



(3) 



where $ is the ionizing luminosity. We note that this hy- 
drodynamic calculation does not include a cold bound 'disc' 
component, since as discussed above the base density struc- 
ture of the wind is known a priori. This base density, along 
with Keplerian rotation is reset at every time-step and the 
model is allowed to evolve to a steady state launching a wind 
from the grid's mid-plane (representing the disc's surface). 
As ex pected, we find e xcell ent agreement wi th the results 
of the lFont etaU l|2004h and I Alexander] (|2008l ) calculations. 
In Figure [T] we show a plot of the converged wind structure 
showing that the wind is approximately spherical once it has 
reached several scale heights, in agreement with our earlier 
discussion in Section [2] 



2.2.1 Numerical EUV Photoevaporative Wind 

In order to determine an accurate kinematic structure of the 
wind we must compute a numerical solution to the prob- 
lem. We use the zeus2d code (Stone & Norman 1992); we 
employ a spherical grid with 6 — [0,7r/2], the radial grid 
cells are logarithmically spaced, such that we have adequate 
resolution at small radius to resolve the onset of the flow. 
The calculation is an isothermal wind calculation with the 
sound speed set to c s — 10km s _1 and the radius scaled 
to the length scale r g = GM*jc\ i.e. the radius at which 
the internal energy of the gas is enough to unbind the gas 
from the star. We use a radial range of r — [0.05, 40]r 9 with 
N r = 240 and Ng = 50. The number density at the base 
(i.e. the density along the 6 — n/2 axis) of the wind was 
calculated semi-analytically by Hollenbach et al. (1994) and 
scales as R~ 3 ^ 2 for R < R g and as R" 5 ^ 2 outside R g , as in 
Font et al. (2004) and Alexander et al. (2 008) we adopt the 
smooth base density profile suggested bv lFont et alj (|2004l ) 
that varies between the two power laws: 

2 1 1/5 

n{R) = 719 [{R/R g ysn + {R/R g ysn\ (2) 



2.3 Dust Distribution 

In order to calculate the dust distribution in the wind, we 
must make some assumptions about the underlying dust dis- 
tribution in the disc. We adopt a dust to gas mass ratio of 
0.01, and a power law grain size distribution of index -3.5 
(Mathis, Rumpl & Nordsiek, 1977, MRN), with grain sizes 
ranging from a m i„=5x 10 _3 /im to a ma x=l mm i we assume 
spherical grains with a density of lg cm -3 . We also assume 
that the dust is fully mixed within the disc up to the transi- 
tion between the bound cold disc and the hot EUV heated 
flow. We then calculate streamlines from the base of the flow 
to the edge of the grid, then along each streamline, compute 
the force balance between the drag force (calculated from 
Equation 1), gravity and the centrifugal force. We take the 
dust grain as entrained if the net force along the streamline 
is > 0. We then obtain the maximum grain size entrained 
along the entire streamline (making sure to check that each 
grain can be entrained the entire length of the streamline). 
In Figure [2] we show the obtained maximum grain size as a 
function of position in the flow for an ionizing luminosity of 
$ = 10 43 S -\ 

We can then compute the dust density from this maxi- 
mum grain size under the assumption that the dust all comes 
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Figure 2. Colour map of the maximum size of entrained dust 
particles at a given position for & = 10 43 s _1 , we note that most 
dust comes from within l-5r g . 
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Figure 3. Dust mass density distribution in the wind model for 
the $ = 10 43 s _1 calculation. 

from the same underlying population described above. In 
Figure [3] we show the obtained dust density, noting that the 
combination of a photoevaporative wind and a selective dust 
population naturally reproduces a 'wingnut' structure. 

2.4 Radiative Transfer 

We calculate the radiative transfer as a two step process: 
First we calculate the temperature structure of the dust us- 
ing a 3D grid calculated from the 2D dust density calcula- 
tions. We have modified the MOCASSIN code (Ercolano et al. 
2003, 2005, 2008) to allow cell-to-cell variations of the grain 
size distribution. We interpolate the results from Section HOI 
onto at 3D Cartesian grid with resolution Nx x Ny X Nz = 
99 x 99 x 50 assuming azimuthal symmetry, using the range 
in the x and y axis of [-32,32]r 3 and in the z axis of [0,32]r 9 . 

The input SED used to irradiate the wind is taken 
from the model set of Robitaille et al. (2006), Model Num- 
ber:3011099, corresponding to a 2.51M0 star with an effec- 
tive temperature of 9480K and radius of 2. 18R0 , surrounded 
by an optically thick disc extending into the dust destruc- 



tion radius. The irradiating SED is placed at the origin of the 
Cartesian grid; this approximation is valid in this case where 
we are mainly interested in scattering of near-IR (NIR) pho- 
tons on the much larger scale of the wind i.e. on scale much 
greater than the scale the disc emission arises from. 

In order to calculate images of the dust distribution, we 
have written a simple ray tracing driver for MOCASSIN, which 
we have called mocassinThinImage. The standard Monte- 
Carlo radiative transfer techniques used to calculate scat- 
tered light images in MOCASSIN (i.e. by capturing scattered 
photons along the line of sight of interest) is inappropriate 
for use in this case. The wind structure is too optically thin 
to obtain a converged image with a computationally feasible 
number of photon packets, since most packets leave the grid 
without interacting with the distribution^. mocassinThin- 
Image performs a ray tracing calculation throughout the 
grid in both scattered and thermal emission, under the as- 
sumption that the medium is optically thin (true of the wind 
always). We note again that since we have not explicitly in- 
cluded the dust in the 'bound' region of the disc we will not 
calculate scattered light emission from this region. We also 
assume the disc mid-plane is optically thick to all radiation 
and hence only consider one hemisphere of the wind. 



2.5 Calculations of Crystallinity Structure 

In order to consider the structure of the dust grains in the 
wind, we have built a very simple dust crystallinity model. 
We assume that the grains are launched from a disc that is 
infinitesimally thin and has the crystalline structure calcu- 
lated by Dullemond et al. (2006), for a Herbig star with a 
mass of 2.5Mq after 3.3Myr of evolution. Any dust that was 
amorphous and is heated to a temperature > 800K while en- 
trained in the wind is then assumed to immediately become 
totally crystalline. Any dust that starts off crystalline or be- 
comes crystalline at any point in the flow is assumed to be 
crystalline for the remainder of the flow regardless of the 
temperature it reaches at any other stage (we find that no 
dust in the flow reaches a temperature exceeding the evap- 
oration temperature, which would destroy crystalline dust 
grains). In order to calculate a radial profile of the fraction 
of dust in the wind that is crystalline/amorphous we must 
apply some cut off on the flow since the total mass of the 
flow is a function of time since it switched on, rather than a 
converged quantity (it diverges logarithmically). Since much 
of the information on dust structure comes from the 10/im 
silicate feature a temperature cut-off is most appropriate, 
thus we calculate the mass-fraction of crystalline grains out 
to a temperature of 100K; as such this should provide a 
reasonable estimate of the observable crystallinity-fraction 
(anything colder is extremely unlikely to make much im- 
pact on the 10/xm silicate feature). Furthermore, since each 
cylindrical radial region will be cut by many streamlines 
with different crystallinity fractions, in order to calculate an 



3 It is important to emphases this is not the case for the tem- 
perature calculation discussed above, where an interaction is not 
required with a photon packet to include its contribution to the 
radiative energy density and hence temperature in that cell (see 
Lucy 1999; Ercolano et al 2003). 
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Table 1. Model Luminosities 
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average crystallinity fraction we perform a number density 
weighted average out to the temperature cut-off. 



3 IMAGES 

In Figure [4] we show calculated images of the 5 different 
models in the H, K' & L' bands, they all show a 'wingnut' 
morphology. The emission is strongest in the bluer bands 
with the H band being strongest and the L' being the weak- 
est. 

As the mass-loss rate increases the emission becomes 
more extended and reaches larger radii, since the density 
along a streamline falls to a given value at larger radius. We 
also find that all of the emission is scattered light from the 
central star and inner disc, many orders of magnitude above 
any thermal emission from the wind itself, in agreement with 
our assertion that the wind is optically thin. The fact that 
the emission is dominated by scattered light allows us to un- 
derstand the evolution of the color variation and emission 
morphology as the mass-loss rates change. For high mass- 
loss rates the emission morphology follows the gas density 
distribution of the wind (i.e. essentially spherical with a gap 
on axis in a region where the wind is negligible) , while more 
pronounced 'wingnut' profiles are clearly visible for lower 
values of mass-loss rate. This arises from the fact that scat- 
tering at a given wavelength is dominated by dust grains of 
comparable sizes. Given the bands we are interested in are 
in the range 1-5/xm then when there is large spatial varia- 
tion in grains of this size then a profile that differs markedly 
from the gas density would be expected (e.g. the wingnut). 
However, for large enough mass-loss rates dust grains with 
sizes in the 1-5/xm range can be entrained in the flow every- 
where and the emission then simply follows the gas density 
distribution. We find that all of the synthetic images are 
dominated by blue light; this is somewhat expected since 
the smallest grains can always be entrained in the flow and 
outnumber the larger grains, thus dominating the opacity 
for an MRN size distribution as the one considered here. 
Furthermore, the negative slope with wavelength of the ir- 
radiating field, which in the 1-5/im range is dominated by 
emission from the disc, means that short wavelength photons 
are more abundant, contributing to the blue appearance of 
the synthetic images. 



4 CRYSTALLINITY PROFILE 

We have used the wind profiles and dust temperature infor- 
mation from the radiative transfer calculation (Section 12. 4p 



to calculate crystallinity fractions in the wind. We have as- 
sessed the impact of the wind on the observable crystallinity 
fraction for discs, where the inner disc is both obscured 
(edge-on) and unobscured. In the unobscured case, to ob- 
serve the disc's crystallinity fraction the wind will also make 
a contribution along the line of sight. We compare the disc's 
emission to the wind's emission at 10/xm and find that at 
all radii the disc dominates over the wind, meaning that the 
observed crystallinity fraction for non-edge on discs is un- 
affected by the photoevaporative wind. However, when the 
disc is observed close to edge-on the disc's contribution will 
be undetectable and the crystallinity fraction in the wind 
maybe studied. 

It is important to note at this point, that the radial 
distribution of the crystallinity fraction in a wind is likely 
to differ from the disc's 'native' distribution. In Figure [5] 
we show radial profiles of the crystallinity fraction obtained 
from a photoevaporative wind (solid black line) and compare 
it to the input crystallinity distribution of the disc (dashed 
blue line). The figure shows that at small radii < 2.5AU 
there is no crystallinity in the wind simply because there is 
no wind inside this radius to carry the dust grains to large 
heights. Between 2.5-200AU there is a falling crystallinity 
fraction; however, it is always in excess of the disc's under- 
lying fraction, whereas outside 200AU it continues to fall 
even though the crystallinity fraction in the disc has now 
risen above the level in the wind. The crystallinity fraction 
in the wind remains the same for all ionizing luminosities 
and is insensitive to the temperature cut-off described in 
Section [2.51 in the range 150-50K. 

The observed wind crystallinity profile is a direct con- 
sequence of the streamline topology, which implies that the 
grains in the wind at any given cylindrical radius must have 
originated in the disc at a smaller radius and hence, in gen- 
eral, with a higher crystallinity fraction. The rise in the 
disc crystallinity at a radii > 100AU is not translated into 
the wind population, since most of the dust in the wind 
is entrained inside 5r 9 « 100AU. The low wind rates at 
> 100AU, coupled with the geometric dilution of the wind, 
mean that the crystallinity fraction in the wind continues 
to fall with radius. The fact that the crystallinity fraction 
is insensitive to the mass-loss rate is simply a consequence 
that the number of dust particles at a given point scales 
approximately as ndust oc $ 1//2 everywhere, the same as the 
gas, since it is largely unaffected by the change in maximum 
grain size for the MRN dust distribution chosen. 



5 COMPARISONS TO OBSERVED EDGE-ON 
OBJECTS 

There are several objects that show extended emission above 
and below the dark dust lane indicative of an optically thick 
disc. All the objects in the sample presented by Padgett et 
al. (1999), which includes the well known 'Butterfly Star', 
show extended scattered light emission above and below a 
dark dust lane. The objects in the Padgett et al. (1999) sam- 
ple; however, all are relatively low-mass young stars and the 
high fluxes detected in scattered light indicate that there 
is much more mass present than would be expected in a 
photoevaporative wind. Furthermore the observed structure 
is much more filamentary than our models predict. These 
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Figure 4. Synthetic images for disc models with irradiating fluxes $ = 10 41 , 10 42 , 10 43 , 10 44 , 10 45 erg s — 1 , the far left hand column 
shows the image in the H band, the next column displays the K band and the next column the L band. The far right column displays 
an RGB composite image (L,K & H bands respectively). The images are individually scaled so that there is a 5dex spread between the 
brightest pixel and the darkest. All images assume the disc is edge-on (i.e. an inclination of 90°), therefore we block out the star assuming 
it would be blocked out by the presence of an optically thick disc. 
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Figure 5. Radial profiles of the wind crystallinity fraction (solid 
line) calculated as a number density weighted mean in the wind 
and the underlying disc distribution taken from Dullcmond et al. 
(2006) for an evolved 3.3Myr 2.5M Herbig star (dashed line). 



observations clearly indicate that a photoevaporative wind 
is too weak to have produced the observed extended emis- 
sion. Padgett et al. (1999) suggested that these objects are 
young (Class I) objects, with the extended emission arising 



from an infalling envelope. Wolf et al. (2003) and Stark et 
al. (2006) have successfully reproduced these observations 
using detailed radiative transfer modelling that include the 
combination of a disc component, an in-falling envelope and 
an outflow cavity. 

PDS 144N (Perrin et al. 2006); however, is a more 
promising candidate for a photoevaporative wind, with im- 
ages showing a morphology very similar to the 'wingnut' 
shapes predicted for the lower mass-loss rates <!> ^ 10 43 erg 
s _1 . Furthermore adopting the <I? = 10 43 erg s _1 , which pre- 
dicts a surface brightness of 0.01 Jy arcsec -2 in H at a radius 
and height of (R = 100, z = 100)AU, we can derive a dis- 
tance to PDS 144N. By comparing the observed image to 
our predicted image we find a distance of 200-400pc, which 
is in the range of previous estimates^. As discussed in Perrin 
et al. (2006) , it is highly unlikely that the morphology arises 
from foreground extinction as suggested explanation to the 
morphology of the Padgett et al. (1999) sample. Perrin et al. 
(2006) indeed discuss photoevaporation as a possible source 
for the extended emission. However, they oversimplify the 
photoevaporation model by just considering the radial scale 
r g , ignoring the fact that, as discussed in Section 2, dust 
grains can be carried to very large distances in the wind, and 



4 Previous estimates range from 140 to 2000pc. Perrin et al. 
(2006) use a nearby A5 Herbig star (PDS 144S) to derive a dis- 
tance of lOOOpc, taking the system to be a wide binary. However, 
such a distance places the system 300pc above the Galactic plane. 
This would require the ejection of the wide binary system from a 
dense cluster, or formation in situ, both of which are unlikely. 
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that at several flow scale heights from the disc the stream- 
lines are spherical. Our models show that the use of r g in 
discussing the characteristic scales of the flow in dust is a 
poor approximation, since, as shown in Figure [2] dust en- 
trained at ~ lr 9 can be easily carried to radii and heights of 
> 20r s producing a dust density that varies on scales differ- 
ent to r g as shown in Figure 3. Perrin et al. (2006) present 
models of the object in which they construct an in-falling 
envelope with a jet cavity surrounding an extended passive 
accretion disc. While their model (Figure 6 of Perrin et al. 
2006) can also reproduce extended emission, it cannot re- 
produce the distinctive 'wingnut' morphology seen in the 
observations and in the model presented in this work (see 
the top panels of Figure [4| . 

However both the model presented by Perrin et al. 
(2006) and the ones presented in this work fail to repro- 
duce one important aspect of the observations. The observed 
colour of PDS 144N is such that the extended regions of the 
emission become dominated by the redder scattered light. 
This is opposite to what is predicted by the Perrin et al. 
(2006) model which predicts a colour variation of red to 
blue with height and the blue band dominating the emis- 
sion at large height. This is due to a dust population whose 
maximum grains size decreases with height, as you move 
from the disc population to the envelope (ISM like) popula- 
tion. On the contrary, our model can reproduce the sign of 
the observed colour variation i.e that the relative strength 
of the red light increases with height. This is demonstrated 
in Figure [6] where we plot the change in H-K' and H-L' for 
the $ = 10 43 s" 1 model as a function of height above the 
mid-plane for the $ = 10 43 s _1 , where there is variations in 
the grain population at sizes of a few microns, which dom- 
inate the scattering opacity at NIR wavelengths. While our 
models predict that the colours become redder with height 
above the disc (due to the fact that larger grains are en- 
trained at greater height, see Figure [5] and these scatter 
red light more efficiently), our models predict scattered blue 
light to remain dominant at large heights (see Figure 4). 
As discussed in Section^ in our model blue light dominates 
due to the spectral slope of the irradiating SED and the fact 
that small grains are present in the wind at all heights. We 
estimate using Figure 3 of Perrin et al. (2006) a change in 
H-L' (A(H-L)) of >1.5, at least three times larger than our 
simple model can predict. However, in order to allow the 
redder bands to dominate over the bluer bands one needs to 
(i) increase the available disc emission at the redder wave- 
lengths so there is more original flux to scatter and/or (ii) 
remove the smaller grains which are currently dominating 
the scattered light (iii) increase the optical depth to the 
star/disc so the extended region becomes optically thick, 
resulting in the reddening of photons as they pass through. 
Grain growth in the disc provides a natural solution to both 
(i) and (ii). Dullemond & Dominik (2005) have shown that 
the dust grains in a disc easily grow to ~ 1/im reducing the 
opacity of the disc at shorter wavelengths and resulting in a 
disc spectrum that falls less steeply or even rises to longer 
wavelengths in the 1-5/xm range. Furthermore, the removal 
of small grains from the disc will also result in the removal 
of small grains from the wind, reducing the scattering effi- 
ciency at the shorter bluer radiation. 

Future observations of the 10/zm silicate feature may 
help disentangle the origin of the dust particles in the ex- 



0.7 



0.6 - 




Height [AU] 

Figure 6. Change in colour (AH-X) as a function of height above 
the midplane, at a cylindrical radius of 100AU. Green line: X=K' 
band colour; red line: X=L' band colour. The model shown is the 
$ = 10 43 s" 1 calculation. 

tended region around PDS 144N since an in-falling envelope 
would be entirely composed of amorphous ISM type grains, 
while, as shown by our models, a photoevaporative wind ori- 
gin implies the presence of crystalline silicate grains in the 
extended regions. 



6 CONCLUSIONS 

We have considered the fate of dust grains that can be en- 
trained in a photoevaporative wind from the surface of a 
disc surrounding a Herbig Ae/Be star. For a median mass- 
loss rate of ~ lO _9 M0yr _1 , we find that grains up to radii 
of several microns can be entrained in the wind. We also 
show that, once entrained in the wind, the dust grains will 
remain entrained and be carried out to very large radius. 

We have considered the observational imprint of this 
wind-entrained dust on edge-on discs, showing that the com- 
bination of a photoevaporative wind structure and a variable 
dust grain population resulting from the variable drag force 
at the base of the wind, can naturally reproduce a 'wingnut' 
morphology of the dust density distribution in the wind. Us- 
ing a combination of Monte-Carlo and ray tracing radiative 
transfer techniques, we calculate scattered light images from 
the density distribution at near infared wavelengths. We find 
our model is not applicable to the well known edge-on discs 
in the Padgett et al. (1999) sample, which are too young 
and optically thick to be explained with a photoevaporative 
origin, but are likely to arise from an infalling envelope as 
shown in Wolf et al. (2003) and Stark et al. (2006). These 
synthetic images show a wingnut morphology inferred from 
the dust density distributions, similar to observations of the 
edge-on disc around PDS 144N (Perrin et al. 2006). The syn- 
thetic images, however, are dominated by emission from the 
smallest grains entrained in the flow, hence failing to repro- 
duce the colour gradient of the observations, which show red- 
der emission at larger heights above the disc. Grain growth 
in the disc, is a natural solution to the colour problem, and 
we estimate severely depleted abundance of grains with radii 
smaller than 1/im. 
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Finally we consider the crystallinity of dust grains en- 
trained in the flow: By following the thermal evolution of 
the grains, we find that crystallinity fraction will remain un- 
changed when the disc's photosphere can be observed since 
the observable mass in the wind in much less than the ob- 
servable mass in the disc. However when the disc photo- 
sphere is obscured (i.e. for edge-on discs) the crystallinity 
fraction in the wind is significantly enhanced above the discs 
underlying crystallinity fraction. Finally, we suggest that de- 
tection of crystalline grains in extended emission around an 
edge on disc is indicative of a photoevaporative wind and 
argues against an envelope origin for the extended emission. 
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